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ABSTRACT 

Interlaminar  damage  in  the  form  of  delamination  due  to  the  stress  concentration  is  one  of  the 
dominant  forms  of  damage  in  laminated  composite  structures.  One  of  the  approaches  to  enhance 
the  delamination  resistance  of  the  advanced  composite  structures  is  to  reduce  mismatch  of  elastic 
properties  and  stress  concentrations  at  the  interfaces  between  the  laminated  layers  by  utilizing 
interleave  materials.  The  interleave  material  considered  in  this  study  is  a  nano-modified  epoxy 
reinforced  with  vapor-grown  nanofibers.  The  nanofibers  are  functionalized  with  ozone  to 
enhance  the  interfacial  strength  between  the  nanofibers  and  the  surrounding  epoxy. 

To  numerically  assess  the  effectiveness  of  this  approach,  continuous  interface  elements  of  zero 
thickness  are  utilized  for  the  characterization  of  the  nano-modified  layer.  The  interface  element 
formulation  is  based  on  a  cohesive  zone  model  with  softening  constitutive  law.  The  numerical 
study  is  carried  out  with  a  focus  on  the  composite  laminates  with  a  stacking  sequence  of  [30/ 
-3 0/90/90/-3 0/3 0]t  under  uniaxial  tensile  loading.  The  stress  analysis  shows  that  the  laminates 
with  such  a  layup  can  have  large  interlaminar  stress  components.  Three-dimensional  finite 
element  analyses  are  carried  out  with  the  interface  elements  being  placed  along  the  ply 
interfaces,  and  the  results  are  compared  with  the  available  experimental  data. 

Computational  and  experimental  studies  were  carried  out  for  the  laminated  composites 
reinforced  with  nano-interlayers  that  were  fabricated  in  the  form  of  either  a  thin  film  of  nano- 
modified  epoxy  reinforced  with  vapor-grown  nanofibers  or  a  fiber  mat  that  can  be  fabricated 
with  the  electrospun  fiber  process.  The  computation  model  used  a  cohesive  zone  model  with 
interface  elements  based  on  bilinear  and  exponential  constitutive  laws.  The  model  was  used  to 
simulate  mode  I,  mode  II  and  mixed-mode  fracture  tests  to  assess  initiation  and  progression  of 
delamination  characteristics  of  the  laminated  composites. 

From  a  numerical  parametric  study  for  the  mode  I  fracture  with  a  double  cantilever  beam  (DCB) 
test,  we  found  that  the  maximum  load  is  mostly  dependent  on  the  critical  energy  release  rate 
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( GIC )  with  the  bilinear  constitutive  law.  Interfacial  strength  and  penalty  parameter  play  fewer 

roles  in  determining  the  maximum  load.  Parametric  studies  were  carried  out  for  mode  II,  an  end 
notch  flexure  (ENF)  test,  to  assess  the  delamination  characteristics  corresponding  to  different 
parameters  of  the  exponential  constitutive  law.  Finally,  uniaxial  tension  simulations  were  carried 
out  for  the  [-30/30/90/NL]s  specimen  to  predict  the  free-edge  delamination  characteristics 
corresponding  to  different  interface  properties. 

Experimental  observation  from  micrographs  and  stress-strain  behavior  shows  that  the  laminated 
composites  with  the  nano-interlayers  can  suppress  the  stress  concentration  and  severity  of  the 
microcrack  and  delamination  damage. 

KEY  WORDS:  Microcracking,  Thermomechanical  Loadings,  Finite  Element  Analysis. 


1  INTRODUCTION 

Interlaminar  damage  in  the  form  of  delamination  due  to  the  stress  concentration  is  one  of  the 
dominant  forms  of  damage  in  laminated  composite  structures.  One  of  the  approaches  to  enhance 
the  delamination  resistance  of  the  advanced  composite  structures  is  to  reduce  mismatch  of  elastic 
properties  and  the  stress  concentrations  at  the  interfaces  between  the  laminated  layers  by 
utilizing  interleave  materials.  The  interleave  material  considered  in  this  study  is  a  nano¬ 
interlayer  which  can  be  fabricated  in  the  form  of  either  a  thin  film  of  nano-modified  epoxy 
reinforced  with  vapor-grown  carbon  nano  fibers  (YGCNFs)  or  a  fiber  mat  that  can  be  fabricated 
with  the  electrospun  fiber  process.  The  nanofibers  in  the  form  of  a  thin  film  are  functionalized 
with  ozone  to  enhance  the  interfacial  strength  between  the  nanofibers  and  the  surrounding  epoxy. 

To  computationally  assess  the  effectiveness  of  this  approach,  continuous  interface  elements  of 
zero  thickness  are  utilized  for  the  characterization  of  the  nano-interlayer.  The  interface  element 
formulation  is  based  on  a  cohesive  zone  model  with  softening  constitutive  law.  Meanwhile,  the 
experiment  was  performed  with  two  types  of  nano-interlayers  that  are  inserted  into  the  laminated 
composites  to  monitor  initiation  and  progression  of  their  damage  behavior.  The  first  section 
illustrates  the  formulation  of  the  computational  cohesive  zone  model  and  the  results  of  the 
parametric  study.  The  next  section  illustrates  the  experimental  observation. 


2  COMPUTATION 

A  cohesive  zone  model  (CZM)  has  gained  much  attention  in  predicting  the  delamination  in 
laminated  composites  structures.  The  CZM  uses  an  interface  element  based  on  a  softening 
constitutive  law.  The  advantage  of  the  CZM  approach  originates  from  the  constitutive  model  of 
the  interface  zone  that  directly  relates  the  opening  displacement  jump  to  the  traction  force  and 
thus  avoids  stress  singularity  problems.  The  CZM  for  predicting  the  initiation  and  progression  of 
the  delamination  in  laminated  composites  was  developed  with  an  interface  element  that  is 
implemented  in  a  finite  element  (FE)  program.  The  model  uses  bilinear  and  exponential 
softening  constitutive  laws  that  describe  the  load  vs.  crack  opening  displacement  behaviors. 
Parallel  implementation  for  solving  nonlinear  finite  element  equations  with  interface  elements 
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was  tested  and  validated.  The  parallelization  of  the  finite  element  program  is  beneficial  both  in 
terms  of  computational  efficiency  and  ability  to  handle  large  problems.  Also,  we  implemented 
and  tested  the  exponential  constitutive  law  which  can  be  more  reliable  to  obtain  the  converged 
solution  than  the  bilinear  constitutive  law.  This  section  illustrates  the  derivation  of  the 
computational  model  and  the  numerical  results  as  well  as  the  parametric  study  with  various  input 
parameters. 

2.1  Constitutive  relations  of  the  interfacial  zone 

The  CZM  incorporates  the  softening  constitutive  behavior  that  relates  the  displacement  jump  to 
the  traction  force  in  the  interface.  The  shape  of  the  softening  constitutive  relation  of  CZM  is 
derived  from  empirical  observations  and  allows  it  to  represent  material  degradation  near  a  crack 
tip  without  having  to  directly  deal  with  the  singular  stress  fields  of  linear  fracture  mechanics. 
Among  several  constitutive  models  adopted  in  developing  the  interface  elements  to  simulate 
delamination,  bilinear  (developed  by  Camanho,  et  al.  [1])  and  exponential  (developed  by  Goyal, 
et  al.  [2])  constitutive  models  are  adopted  in  the  present  study.  Both  models  are  capable  of 
handling  mixed-mode  delamination  and  account  for  the  irreversibility  of  the  damage  process. 
Equations  for  the  two  constitutive  models  are  given  in  the  remainder  of  this  section. 

To  achieve  accurate  prediction  of  onset  and  growth  of  the  delamination,  especially  for  complex 
three-dimensional  composite  structures,  sufficiently  refined  FE  meshes  are  required  especially  in 
the  direction  of  delamination  propagation.  Furthermore,  the  class  of  nonlinear  problems 
involving  softening  constitutive  law  often  exhibit  sharp  nonlinearities  and  requires  many 
iterations  with  the  nonlinear  solution  methods  such  as  a  Newton-Raphson  method  to  obtain 
convergence  at  each  load  step.  Therefore,  a  high-performance  computing  technology  is 
harnessed  to  carry  out  the  computer-intensive  nonlinear  analysis  efficiently.  The  interface 
elements  are  implemented  in  a  parallel  FE  program.  It  is  characterized  by  incorporation  of  a 
memory-  and  computing-efficient  multifrontal  solver  [3]  in  an  object-oriented  interpreter  layer  in 
a  Python  computer  programming  language  and  parallel  computation  of  residual  force  vectors 
during  iterations. 


2.2  Numerical  results  and  parametric  study 
Mode  I  f  racture 

We  predicted  a  model  I  fracture  behavior  with  the  CZM  model.  The  model  simulated  a  double 
cantilever  beam  (DCB)  test  configuration  as  shown  in  Figure  1.  The  stacking  sequence  of  the 
lamination  is  [+30/-30/90/NL/90/-30/+30]T  where  NL  is  a  nano-interlayer  with  prescribed 
partially  delaminated  zone  between  the  [90/90]  interface  surface.  The  stress  analysis  shows  that 
the  laminates  with  such  a  layup  can  have  large  interlaminar  stress  components  that  contribute  to 
the  delamination  damage.  The  composite  layer  is  made  with  a  carbon-fiber  reinforced  epoxy 
material  whose  properties  are  listed  in  Table  1. 
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Figure  1. 

Configuration  of  a 
double-cantilever 
beam  test. 


Table  1 .  Material  properties  of  a  carbon- fiber  reinforced  epoxy  material 


Eu 

^22  —  ^33 

Gu=Gn 

G23 

n2=n3 

179.95 

GPa 

9.856 

GPa 

6.0674 

GPa 

3.3998 

GPa 

0.37 

V23 

Gic 

GIIC 

T3 

o 

Tl  ~T2 

0.45 

- 

~ 

The  first  case  considered  is  the  bilinear  softening  constitutive  law.  Several  parameters  determine 
the  shape  of  the  bilinear  law,  and  these  are  interlaminar  strength  ( r3 ),  penalty  parameter  ( K ), 
critical  strain  energy  release  rate  ( GIC ),  critical  stretching  length  (Ac)  and  final  stretching  length 
( Ap  ).  z'l  and  Ac  are  the  load  and  displacement  at  the  peak  point  of  the  bilinear  curve, 
respectively,  if  is  a  slope  of  the  increasing  portion  of  the  line.  GIC  is  the  area  under  the 
bilinear  curve.  r3  and  Ac  are  related  with  the  initiation  of  the  delamination,  and  GIC  is  related 
with  progression  of  the  delamination.  Note  that  a  ratio  Ac/Ap  indicates  brittleness  of  the 

material.  Among  these  parameters,  not  all  of  them  are  independent.  For  example,  the  penalty 
parameter  is  related  with  the  interlaminar  strength  and  the  critical  stretching  length  by 
K  =  t°  / Ac ,  and  the  strain  energy  release  rate  is  expressed  with  GIC  =  GAp/2. 

Three  different  bilinear  softening  constitutive  laws  were  used  for  the  simulation  for  a  parametric 
study.  The  three  cases  are 

I.  same  K  and  Ap ,  different  r" ,  Ac  and  GIC  (see  Figure  2  (a)), 

II.  same  K  and  GIC ,  different  zf ,  Ac  and  Ap  (see  Figure  2  (b)) 

III.  same  Ac  and  A p ,  different  K ,  r3  and  GIC  (see  Figure  2  (c)). 

Figure  2  shows  the  bilinear  laws  of  three  cases.  On  the  right  side  of  the  figures  are  tables  that  list 
three  different  values  of  the  parameters.  Note  that  cases  I  and  II  use  K  =  106  N/mm3. 


4 


GIC  (J/m2) 

t,°  [MPa] 

(a) 

183.8 

28 

(b) 

262.6 

40 

(c) 

341.4 

52 

GIC  (J/m2) 

V  [MPa] 

(a) 

262.6 

28 

(b) 

40 

(c) 

52 

G|C  (J/m2) 

V  [MPa] 

^[N/mm3] 

(a) 

183.8 

28 

0.7x1 06 

(b) 

262.6 

40 

I.OxlO6 

(c) 

341.4 

52 

1.3x106 

Figure  2.  Three  different  bilinear  softening  constitutive  laws  with  three  different  cases,  (a)  Case 
I  with  same  K  and  Ap ,  different  r3 ,  Ac  and  GIC ,  (b)  case  II  with  same  K  and  G1C , 

different  r3 ,  Ac  and  Ap  and  (c)  case  III  with  same  Ac  and  Ap ,  different  K ,  r°  and 


Figure  3  shows  load-displacement  curves  under  the  DCB  tests  by  using  the  three  different 
bilinear  constitutive  laws  in  Figure  2.  The  maximum  peak  loads  ( P )  in  the  figures  are 
dependent  on  the  parameters  varied  in  the  parametric  study.  From  Figure  3  (a)  and  (b),  higher 
r3  and  GIC  results  in  higher  maximum  loads.  It  is  interesting  to  note  that  case  (2),  which  is 

under  the  condition  of  the  same  K  and  GIC  but  different  r3  ,  Ac  and  Ap ,  results  in  nearly 
identical  maximum  loads.  Therefore,  we  can  conclude  that  the  maximum  load  is  not  sensitive  to 
r3  ,  Ac  and  A  p  but  highly  dependent  on  K  and/or  GIC .  The  independence  of  the  maximum 

load  on  r3  and  Ac  indicates  that  the  global  behavior  of  the  DCB  specimen  under  the  mode  I 
fracture  loading  is  not  sensitive  to  the  initiation  of  the  delamination  damage. 
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(a) 


(b) 


(c) 


Figure  3.  Prediction  of  mode  I  fracture  under  a  DCB  test  with  three  different  cases,  (a)  Case  I  in 
Figure  2  (a),  (b)  case  II  in  Figure  2  (b)  and  (c)  case  III  in  Figure  2  (c). 
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Figure  4  shows  a  comparison  of  the  results  from  case  I  and  case  III,  which  describes  the  bilinear 
constitutive  laws  with  the  same  initial  slope  ( K )  and  the  same  starting  stretching  displacement 
( Ac ),  respectively.  Both  the  maximum  load  vs.  GIC  in  Figure  4  (a)  and  the  maximum  load  vs. 

r3°  in  Figure  4  (b)  curves  of  cases  I  and  III  are  identical.  Therefore,  we  can  conclude  that  the 
maximum  load  is  not  sensitive  to  K ,  Ac  or  the  shape  of  the  bilinear  curves.  From  this 
conclusion  and  the  one  observed  earlier,  we  can  conclude  that  the  maximum  load  is  mostly 
dependent  on  the  critical  strain  energy  release  rate  ( GIC ).  The  interfacial  strength  and  penalty 
parameter  play  fewer  roles  in  determining  the  maximum  load.  This  is  an  important  observation 
because  we  only  need  to  pay  attention  to  obtain  the  accurate  GIC  from  the  experiment,  especially 
when  measuring  the  interfacial  strength  is  a  tough  task. 


Figure  4.  Comparison  of  case  I  and  case  III  in  Figure  2.  (a)  max.  load  vs.  GIC  and  (b)  max.  load 
vs.  t'l .  Same  K  and  same  Ac  indicate  case  I  and  III,  respectively. 


Mode  U  fracture 

For  mode  II  delamination,  we  simulated  an  ENF  model  using  the  exponential  cohesive 
constitutive  law  (see  Figure  5).  An  isotropic  material  property  with  E  =  150  GPa  and  v  =0.25 
was  used  for  the  numerical  analyses.  The  ENF  specimen  was  discretized  using  100  elements  in 
the  axial  direction,  and  the  interface  elements  with  exponential  CZM  were  placed  in  the 
interlayer. 
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1.5  mm  x  2 
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Figure  5. 

Configuration  of  an 
ENF  test. 
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Several  parameters  determine  the  shape  of  exponential  constitutive  law.  In  case  of  mode  II 
failure,  these  include  an  interfacial  strength  ( r , 0 )  which  defines  the  peak  value  of  the  traction- 
displacement  curve,  a  critical  energy  release  rate  ( Gnc )  which  is  the  area  under  the  curve,  and  a 
shape  parameter  (/?)  which  determines  the  general  shape  of  the  exponential  curve.  Figure  6 
shows  the  traction-displacement  curves  with  three  different  shape  parameters.  Parametric  studies 
were  conducted  to  see  the  effect  of  these  parameters  on  the  mode  II  response  of  the  ENF  test 
specimen. 


Figure  6. 

Traction-displacement 
curves  of  exponential 
constitutive  model 
with  three  different 
shape  parameters  (/?). 


Figure  7  shows  the  load-displacement  response  and  the  initial  peak  load  of  the  ENF  specimens 
for  different  values  ranging  of  GIIC  from  0.68  N/mm  to  2.32  N/mm.  Correlation  between  initial 

peak  load  and  Gnc  can  be  observed  for  the  Mode  II  failure  with  the  exponential  model. 
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Figure  7.  Response  of  the  ENF  specimens;  (a)  load  displacement  curve,  and  (b)  initial  peak  load 
for  different  Gnc  values. 


Plots  in  Figure  8  depict  the  effect  of  different  interfacial  strength  and  /?  on  the  initial  peak  load 
value.  The  critical  energy  release  rate  is  kept  constant  ( G,IC  =  1 .25  N/mm )  for  the  parametric 
studies.  While  the  previous  results  of  Mode  I  with  the  bilinear  constitutive  model  showed  peak 
load  was  determined  by  GIC ,  different  peak  values  could  be  observed  by  varying  the  parameters 

that  determine  the  shape  of  the  exponential  constitutive  model  although  a  constant  Guc  was  used 
for  the  ENF  simulation. 


Figure  8.  Initial  peak  load  of  ENF  specimens;  (a)  peak  load  vs.  r,"  and  (b)  peak  load  vs.  /? 
Mixed-mode  fracture 

The  composite  laminate  with  a  stacking  sequence  of  [+30/-30/90/NL/90/-30/+30]t  is  subjected 
to  uniaxial  tension  loading.  The  test  specimen  has  a  length,  width  and  thickness  of  254  mm,  25.4 
mm  and  8.382  mm,  respectively. 

For  the  FE  analysis,  the  specimen  was  discretized  with  120x30x6  8-node  hexahedral  elements 
(90024  DOFs),  and  the  interface  elements  with  the  exponential  CZM  (/?  =  1.5)  were  placed  in 
the  middle  of  the  [90/90]  interlayer  to  simulate  the  initiation  and  propagation  of  the  delamination. 
Since  the  damage  is  mainly  caused  by  the  stress  concentration  in  the  ffee-edge  region,  the 
specimen  was  modeled  in  such  a  way  that  the  FE  mesh  was  refined  from  the  center  toward  the 
free  edges.  The  displacement-controlled  Newton-Raphson  technique  with  automatic  increment 
adjustment  and  the  line  search  algorithm  were  employed  to  obtain  the  converged  solution  of  the 
nonlinear  equations. 


9 


Figure  13  shows  the  evolution  of  delamination  damage  represented  by  the  contour  of  damage 
variable  with  the  strain  increase.  The  value  of  damage  variables  greater  than  zero  and  equal  to 
one  imply  the  initiation  and  propagation  of  the  delamination,  respectively. 


10 


Figure  9.  Contours  of  damage  at  various  strain  levels  of  (a)  e  =  0.55%  ,  (b)  s  =  0.69%  and 
(c)  5  =  1.18%. 


The  load  and  strain  values  at  which  delamination  initiates  (Tc  and  sc)  and  complete  detachment 
(Tp  and  sp)  occurs  for  different  critical  energy  release  rates  are  summarized  in  Table  2. 


Table  2.  Load/strain  of  delamination  initiation/propagation 


G,c  (N/m) 

GIIC  (N/m) 

Tc(N) 

Sc(%) 

TP(N) 

sp(%) 

184 

490 

6200 

0.43 

8640 

0.65 

263 

701 

7042 

0.49 

9800 

0.70 

328 

911 

7900 

0.55 

11700 

0.83 

3  EXPERIMENT 

Experiment  was  performed  with  the  laminated  composites  reinforced  with  nano-interlayers  to 
evaluate  the  microcrack-  and  delamination-resistance  performance.  We  considered  the  cases 
with  the  nano-interlayers  made  in  the  form  of  (a)  thin  film  of  epoxy  material  reinforced  with  the 
VGCNFs  and  (b)  fiber  mat  that  was  fabricated  with  the  electrospun  process. 

3.1  Film  nano-interlayer 

The  baseline  layup  is  [+30/-30/90]s,  where  the  subscript  s  represents  a  symmetric  layup.  With 
the  baseline  layup,  one,  three  and  five  nano-interlayers  can  be  placed  between  the  composite 
layers  to  make  [+30/-30/90/NL]s,  [+30/-30/NL/90/NL]s,  [+30/NL/-30/NL/90/NL]s,  respectively, 
where  NL  represents  the  nano-interlayer.  The  nano-interlayer  was  fabricated  in  the  form  of  thin 
film  with  epoxy  material  reinforced  with  VGCNFs. 

The  laminated  composites  of  the  [+30/-30/90]s  layup  were  loaded  under  uniaxial  tension  loading 
at  a  stress  level  of  50  ksi.  The  free  edge  of  the  test  specimen  was  polished  and  observed  with  a 
microscope.  Figure  10  shows  the  micrographs  taken  at  the  free  edge  of  the  specimen  with  and 
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without  the  nano-interlayer.  As  the  analysis  predicts,  severe  interlaminar  stress  causes 
delamination  in  the  middle  of  the  test  specimen  between  the  [90/90]  layup,  which  thus  causes 
extensive  delamination,  as  shown  in  a  micrograph  in  Figure  10  (a).  With  the  presence  of  the 
nano-interlayer,  the  interlaminar  stress  concentration  is  relieved  to  help  the  reduction  of  the 
severity  of  delamination  damage  as  shown  in  Figure  10  (b). 

Figure  1 1  (a)  shows  the  numbers  of  microcracks  in  a  6-inch  span  at  various  stress  levels  under 
uniaxial  tension  loading.  The  hollow  rectangle  and  solid  circular  markers  represent  the  crack 
counts  without  and  with  the  nano-interlayers,  respectively.  The  figure  clearly  shows  that 
existence  of  the  nano-interlayer  in  the  laminated  composites  suppresses  the  microcracks  with  the 
increase  of  the  loading.  Figure  1 1  (b)  shows  the  stress  level  at  the  initiation  of  the  delamination 
without  and  with  the  nano-interlayers.  We  can  observe  a  significant  increase  of  the  delamination 
initiation  stress  with  the  presence  of  the  nano-interlayer. 


Figure  10.  Micrographs  with  microcracks  and  delamination  under  a  uniaxial  tension  loading  at  a 
stress  level  of  50  ksi:  (a)  no  nano-interlayer  [+30/-30/90]s  and  (b)  with  nano-interlayer 
[+30/-30/NL/90/NL]s. 


(b)  60  t 
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Figure  1 1  (a)  Number  of  microcracks  counted  in  a  6-inch  span  in  the  middle  of  the  test  specimen 
at  various  stress  levels  and  (b)  delamination  stress  with  and  without  the  nano¬ 
interlayers  under  uniaxial  tension  loading. 


3.2  Electrospun  nano-interlayer 

The  present  work  considers  the  electrospun  fiber  mat  as  the  nano-interlayer  to  relieve  the 
interlaminar  stress  concentration.  The  fiber  mats  were  fabricated  with  the  in-house  electrospun 
machine  built  at  Soongsil  University  in  Korea.  The  fiber  mats  were  delivered  to  the  University 
of  Dayton  Research  Institute  (UDRI)  to  be  ready  to  put  between  the  composite  layers.  Once 
fabricated  and  processed,  the  test  specimens  with  layups  of  [+30/(NL)/-30/(NL)/90/(NL)]s  were 
prepared  for  similar  uniaxial  tension  tests  that  were  explained  earlier  in  a  previous  section,  where 
(NL)  indicates  a  possible  location  of  the  insertion  of  the  nano-interlayer.  Therefore,  five  nano¬ 
interlayers  were  inserted  between  the  laminate  plies.  As  shown  in  the  numerical  calculation, 
these  nano-interlayers  are  expected  to  suppress/delay  the  delamination  damage  by  reducing  the 
interlaminar  stress  concentration  at  the  free  edge. 

The  test  specimens  were  loaded  under  uniaxial  tension  with  the  MTS  machine.  An  acoustic 
emission  device  was  used  to  monitor  the  damage  events  due  to  the  microcracks  and  delamination 
during  the  loading.  Once  the  events  were  observed,  we  stopped  the  load  increase  and  monitored 
the  microcracking  and  delamination  damage  by  using  an  optical  microscope.  At  this  time,  we 
counted  the  number  of  microcracks  and  delamination  and  took  some  micrographs.  The 
micrographs  were  taken  for  the  test  specimens  before  loading  (Figure  12),  at  the  first-ply  failure 
(Figure  13)  and  after  the  first  observation  of  delamination  (Figure  14).  In  these  figures,  the 
micrographs  on  the  left-hand-side  and  right-hand  side  represent  the  specimens  without  and  with 
the  nano-interlayers,  respectively.  Figure  12  clearly  shows  the  existence  of  the  nano-interlayers. 
The  thicknesses  measured  from  five  specimens  without  the  nano-interlayers  and  five  specimens 
with  the  five  nano-interlayers  are  approximately  1.229  mm  (0.0484  in)  and  1.229  mm  (0.0484 
in),  respectively.  Therefore,  the  total  thickness  of  the  five  nano-interlayers  from  the  electrospun 
mat  is  less  than  0.001  mm,  which  is  significantly  thinner  than  the  other  existing  interleave  layers. 

Figure  13  shows  the  single  vertical  microcrack  that  is  observed  at  the  free  edge  of  the  specimens 
without  and  with  the  nano-interlayers  after  they  were  loaded  to  165.2  MPa  (24.0  ksi)  and  171.2 
MPa  (24.8  ksi),  respectively.  The  averages  of  the  FPF  stress  with  the  five  pristine  and  five  nano 
specimens  are  175.3  MPa  (25.4  ksi)  and  184.4  MPa  (26.8  ksi),  respectively.  Therefore,  we 
observed  a  5.2%  increase  of  the  FPF  stress  with  the  addition  of  the  nano-interlayers.  Figure  14 
shows  both  the  vertical  microcrack  and  the  horizontal  delamination  that  are  observed  at  the  free 
edge  of  the  specimens  without  and  with  the  nano-interlayers  after  they  were  loaded  to  246.0  MPa 
(35.8  ksi)  and  271.1  MPa  (39.3  ksi),  respectively.  The  averages  of  the  delamination  stress  with 
the  five  pristine  and  five  nano  specimens  are  249.1  MPa  (36.1  ksi)  and  269.1  MPa  (39.0  ksi), 
respectively.  Therefore,  we  observed  an  8.0%  increase  of  the  delamination  stress  with  the 
addition  of  the  nano-interlayers.  The  averages  of  the  numbers  of  microcracks  at  the 
delamination  stress  with  the  five  pristine  and  five  nano  specimens  are  64  and  50,  respectively, 
which  is  a  21.6%  decrease  with  the  addition  of  the  nano-interlayers.  Therefore,  we  can  conclude 
that  the  nano-interlayer  can  suppress  the  microcracking  and  delamination  damage.  The  averages 
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of  the  ultimate  strength  with  the  five  pristine  and  five  nano  specimens  are  296.0  MPa  (42.9  ksi) 
and  324.9  MPa  (47.1  ksi),  respectively.  Therefore,  we  observed  a  9.8%  increase  of  the  ultimate 
strength  with  the  addition  of  the  nano-interlayer. 


Figure  12.  Micrographs  of  test  specimens  before  loading:  (a)  without  nano-interlayer  [+30/ 
-30/90]s  and  (b)  with  nano-interlayer  [+30/NL/-30/NL/90/NL]s. 
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Figure  13.  Micrographs  with  a  microcrack  at  the  first-ply  failure  under  a  uniaxial  tension  loading 
at  a  stress  level  of  (a)  165.2  MPa  (24.0  ksi)  without  nano-interlayer  [+30/-30/90]s  and 
(b)  171.2  MPa  (24.8  ksi)  with  nano-interlayer  [+30/-30/NL/90/NL]s. 
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Figure  14.  Micrographs  after  the  first  observation  of  delamination  under  a  uniaxial  tension 
loading  at  a  stress  level  of  (a)  246.0  MPa  (35.8  ksi)  without  nano-interlayer  [+30/ 
-30/90]s  and  (b)  261.5  MPa  (37.9  ksi)  with  nano-interlayer  [+30/-30/NL/90/NL]s. 


Figure  15(a)  shows  the  number  of  the  microcracks  counted  in  a  6-inch  span  at  the  free  edge  in 
the  middle  of  the  test  specimen  at  various  stress  levels.  The  hollow  dots  with  a  dotted  line  and 
solid  dots  with  a  solid  line  represent  the  data  from  the  five  specimens  without  and  with  the  nano¬ 
interlayer,  respectively.  The  figure  clearly  shows  that  the  number  of  microcracks  is  decreased 
with  the  addition  of  the  nano-interlayers  at  several  stress  levels. 

Figure  15(b)  shows  stress  levels  at  the  first-ply  failure,  delamination  and  ultimate  failure  with 
and  without  the  nano-interlayers  under  uniaxial  tension  loading.  As  stated  above,  the  stress 
levels  at  the  first-ply  failure,  delamination  and  ultimate  failure  increased  5.2%,  8.0%  and  9.8% 
with  the  addition  of  the  nano-interlayers,  respectively.  Therefore,  we  can  conclude  that  the 
presence  of  the  nano-interlayer  can  decrease  and/or  suppress  the  microcracking  and  delamination 
damage  by  reducing  the  interlaminar  stress  concentration  at  the  free  edges  in  the  laminated 
composites. 
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Figure  15.  (a)  Number  of  microcracks  counted  in  a  6-inch  span  at  the  free  edge  in  the  middle  of 
the  test  specimen  at  various  stress  levels  and  (b)  stress  levels  at  the  first-ply  failure, 
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delamination  and  ultimate  failure  with  and  without  the  nano-interlayers  under 
uniaxial  tension  loading. 


4  CONCLUSIONS 

Computational  and  experimental  studies  were  carried  out  for  the  laminated  composites 
reinforced  with  nano-interlayers  that  were  fabricated  in  the  form  of  either  a  thin  fdm  of  nano- 
modified  epoxy  reinforced  with  YGCNFs  or  a  fiber  mat  that  can  be  fabricated  with  the 
electrospun  fiber  process.  The  computation  model  used  a  CZM  with  interface  elements  based  on 
bilinear  and  exponential  constitutive  laws.  The  model  was  used  to  simulate  mode  I,  mode  II  and 
mixed-mode  fracture  tests  to  assess  initiation  and  progression  of  delamination  characteristics  of 
the  laminated  composites. 

From  a  numerical  parametric  study  for  the  mode  I  fracture  with  a  DCB  test,  we  found  that  the 
maximum  load  is  mostly  dependent  on  the  critical  energy  release  rate  ( GIC )  with  the  bilinear 

constitutive  law.  Interfacial  strength  and  penalty  parameter  play  fewer  roles  in  determining  the 
maximum  load.  Parametric  studies  were  carried  out  for  mode  II,  an  ENF  test,  to  assess  the 
delamination  characteristics  corresponding  to  different  parameters  of  the  exponential  constitutive 
law.  Finally,  uniaxial  tension  simulations  were  carried  out  for  the  [-30/30/90/NL]s  specimen  to 
predict  the  free-edge  delamination  characteristics  corresponding  to  different  interface  properties. 

Experimental  observation  from  micrographs  and  stress-strain  behavior  shows  that  the  laminated 
composites  with  nano-interlayers  can  suppress  the  stress  concentration  and  severity  of  the 
microcrack  and  delamination  damage.  In  addition  to  the  potential  enhancement  of  the 
multifunctionality  due  to  its  superior  material  properties,  such  as  thermal  and/or  electrical 
conductivity,  less  thermal  mismatch,  etc.,  the  improvement  of  the  mechanical  performance  with 
the  nano-interlayer  draws  further  investigation  of  the  nano-interlayer  study. 
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